A study of carbon substitutions in MgB2 within the two-band Ehashberg theory 
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We study the effects of C substitutions in MgB2 within the two-band model in the Eliashberg 
formulation. We use as input the B — B stretching-mode frequency and the partial densities of states 
N^{Ef) and N^{Ef), recently calculated for Mg(Bi_xCx)2 at various x values from first-principles 
density functional methods. We then take the prefactor in the Coulomb pseudopotential matrix, fi, 
and the interband scattering parameter, V'^ , as the only adjustable parameters. The dependence on 
the C content of Tc and of the gaps (Ao- and A^) recently measured in Mg(Bi_xCx)2 single crystals 
indicate an almost linear decrease of fi on increasing x, with an increase in interband scattering 
that makes the gaps merge ai x = 0.132. In polycrystals, instead, where the gap merging is not 
observed, no interband scattering is required to fit the experimental data. 

PACS numbers: 74.45.-l-c, 74.70.Ad, 74.20.Fg 



In spite of its simple structure, the intermetallic com- 
pound MgB2 - discovered to be superconducting at about 
40 K in 2001 0] - soon revealed a number of surprising 
features that could not be explained within a picture of 
conventional superconductivity. Bandstructure calcula- 
tions ,;2j showed that the energy bands of MgB2 can be 
grouped into two sets: the quasi-2D a bands, and the 
3D TT bands, originating from the superposition of in- 
plane and out-of-plane boron orbitals, respectively. As 
a matter of fact, most of the physical properties of this 
superconductor have found a clear and relatively simple 
explanation within an effective two-band model |^ ^ |^ 
in which the two bands interact via a phonon-mediated 
interband coupling. The result is that superconductivity 
develops in both bands at the same Tc, but with energy 
gaps of different amplitude, A^- and A^r, and thus dif- 
ferent values of the gap ratio 2A//cb?c. The success of 
the two-band model in describing the features of MgB2 
naturally opens the question whether it can predict (or 
at least explain a posteriori) the effects of induced disor- 
der, irradiation and, over all, chemical substitutions on 
the physical properties of the compound. As far as sub- 
stitutions are concerned, the experimental test of theo- 
retical predictions has been delayed or even prevented by 
the technical difficulties in obtaining good-quality sam- 
ples of partially substituted MgB2 p\ . Recently, point- 
contact measurements of the gap amplitudes as a func- 
tion of the C content have been reported in state-of-the- 
art Mg(Bi_xCx)2 polycrystals and single crystals 
The availability of these results (that for some aspects 
contrast with each other) gives a good opportunity to 
test the two-band model. In this paper we will show 
that both the experimental data concerning Tc and the 
gaps as a function of x can be well explained within the 
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two-band model in the Eliashberg formulation. We will 
use as input the frequencies of the B-B stretching mode 
(which is strongly coupled to the holes in the a band) and 
the partial densities of states at the Fermi level, N^{Ep) 
and N^(E-p), calculated from first-principle density func- 
tional methods adopting the viewpoint of ordered su- 
percells instead of the virtual-crystal approximation. 
Then, we will show that the experimental x dependence 
of Tc and of the gaps A^ and A^ can be very well re- 
produced by admitting a reasonable x dependence of the 
prefactor in the Coulomb pseudopotential matrix 0, Hoj 
and, in the case of single crystals, an increase in the in- 
terband scattering T°"^ on increasing the C content. 

Let us start from the generalization of the Eliashberg 
theory [lllll^ for systems with two bands Ts^, that has 
already been used with success to study the MgB2 system 
H H 113, El m Ell. To obtain the gaps and the criti- 
cal temperature within the s-wave, two-band Eliashberg 
model one has to solve four coupled integral equations 
for the gaps Ai(iLj„) and the renormalization functions 
Zi{iujn), where i is a band index and Wn are the Matsub- 
ara frequencies. We included in the equations (explicitly 
reported elsewhere E3) the non-magnetic impurity scat- 
tering rates in the Born approximation, P^. 

The solution of the Eliashberg equations requires as in- 
put: i) the four (but only three independent jlJl) electron- 
phonon spectral functions afj{uj)F{uj)] ii) the four (but 
only three independent 13] ) elements of the Coulomb 
pseudopotential matrix ^,*{u)c)', iii) the two (but only one 
independent pjj) effective impurity scattering rates F'j. 
None of these parameters or functions has been calcu- 
lated for C-substituted MgB2, and in many cases their 
determination is a very difficult task, at least at the 
present moment. Hence, we are left with a problem with 
too many degrees of freedom. However, we will now show 
how some reasonable approximations allow reducing the 
number of adjustable parameters to 2, with no significant 
loss of generality. 
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FIG. 1: (a) Calculated density of states at the Fermi energy, 
N^{Ef) (solid line) and N^{Ef) (dashed line) as a function 
of X (From Ref. |^. (b) Calculated frequency of the B-B bond 
stretching mode (the i?2g mode in pure MgB2) as a function 
of a; (from Ref. 



Let's start with the four spectral functions a| (a;)i^(ti'), 
that were calculated for pure MgB2 in ref. ll(T For sim- 
plicity, we will assume that the shape of the afjF{uj,x) 
functions does not change with x, and we will only rescale 
them with the electron-phonon coupling constants Ay : 



ayF(w,x) 



A.ij(x) 
Ay(x = 0)' 



,F(w,x = 0) (1) 



Neglecting the effect of C substitution on the shape of the 
e-ph spectral functions is not a dramatic simplification, 
since we showed in a previous paper Q that the details 
of a^F(ijj) do not significantly affect the resulting Tc- To 
determine the scaling factor in eq. ^ let us remind that, 
from the definition of electron-phonon coupling constant 

IF 



m 



A = 



(2) 



where M is the ion mass, J^o is a frequency representa- 
tive of the phonon spectrum, N^{E-p) is the density of 
states at the Fermi level and < > is the average ma- 
trix element of the electron- ion interaction [T^ . In our 
case, M is the boron mass 3] and does not depend on 
X. As a first ap prox imation, and as we did in the case of 
Al substitution [lj| , we will assume that also the average 
matrix element of the electron-ion interaction < > is 
constant, because it is basically determined by the defor- 
mation potential which is almost independent of x |l8j |. 
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FIG. 2: Electron-phonon coupling constants \ij calculated as 
a function of x according to eqs.|3and|l| 



The partial densities of states at the Fermi level, N^{Ep) 
and N^{Ef), have been recently calculated from first 
principles by using a supercell approach for different 
values of the C content x, and are reported in Fig^a). 
The frequency f2o can be identified with the frequency of 
the B-B bond-stretching phonon mode (the E2g mode), 
that has been recently calculated as a function of x from 
first principles 9], and is reported in Fig^b). Since this 
mode couples strongly with the holes on top of the a 
band, from eq. |21 we will have for X^ra (which gives the 
most important contribution to superconductivity in our 
system): 



Act (J (x) 



N-iEF,x)ulJx^O) 



N^iEF.x 



A,,(x = 0). (3) 



In this way, we assume that the change in the frequency 
of the i?2g phonon affects the coupling constant, while we 
neglect its influence on the shape of the electron-phonon 
spectral function. For the other coupling constants, we 
will instead assume for simplicity 



V(j, j) 7^ (o-,cr) Aij(a;) 



N^{Ef,x^O) 



Aij(a: = 0) (4) 



with X^^{x ^ 0)^1.017, X^^{x = 0)^0.448, X„^{x = 
0)=0.2I3 and X^„{x = 0)=0.155 H |l3. Fig. H shows 
the calculated electron-phonon coupling constants A^ as 
a function of x. 

As far as the Coulomb pseudopotential is concerned, 
let us start from its expression in pure MgB2 0, ITol IT^ : 



^,*{x)-^ 



--^l{LOc,x)N^°\EF,x) 



2.23 



1 2.48 
rtot / 



(5) 



where /x(wc, x) is a free parameter and N^^{Ef, x) is the 
total normal density of states at the Fermi level. The 
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FIG. 3: The experimental Tc measured in Mg(Bi_xCx)2 sin- 
gle crystals I2H (solid circles) and polycrystals (open cir- 
cles) as a function of x. The line is only a guide to the eye. 

numbers 2.23 and 2.48 in the Coulomb matrix have been 
calculated for pure MgB2 but, as a first approximation, 
we will suppose them not to depend on x. In this way, 
the elements of the Coulomb pseudopotential matrix, /.t*j, 
depend on x only through the densities of states at the 
Fermi level and through the common prefactor fj,{ujc,x), 
that could also take into account the effects of disorder. 

As far as the scattering rates are concerned, let us re- 
mind that, due to Anderson's theorem, intrahand scatter- 
ing does not affect neither Tc nor the gaps 2(|, so we will 
disregard both F'^'^ and F'^'^. The remaining interband 
scattering parameters are related to each other since 

Aij(x) ^ MyW ^ ^ nI{Ef,x) 

Aji(x) ^Jil{x) Vy^ix) N^iEF,x) 

and thus we will always refer only to F'^'^. Finally, we 
can fix the cut-off energy (e.g., uJc — 700 meV) so as to 
reduce the number of adjustable parameters to two: the 
prefactor in the Coulomb pseudopotential, /i(wc) (that 
we will call simply fi from now on) and the interband 
scattering parameter T"'^ . 

As already pointed out, the aim of the present work is 
to show that the experimental dependence of Tc and of 
the gaps. Act and A^, on the C content in Mg(Bi_xCx)2 
can be explained within the two-band Eliashberg theory. 
The experimental Tc{x) curves measured in single crys- 
tals |2J and polycrystals are reported in Fig|3| The 
corresponding x dependencies of the gaps measured by 
point-contact spectroscopy (PCS) are reported in Fig. ^ 
and Fig. |S1 respectively (symbols). In single crystals 
(Fig. 2)), the two gaps approach each other on increas- 
ing X, until at X = 0.132 they become experimentally 
indistinguishable. This means that, at this doping con- 
tent, their amplitudes are equal to each other within the 
experimental uncertainty. In polycrystals, instead, the 
two gaps remain clearly distinct up to x = 0.10, where 
A^ is much smaller than in single crystals with the same 
C content (see Fig. (S)). 
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FIG. 4: Symbols: experimental values of Aa (open circles) 
and A,r (solid circles) measured at T = 4.2 K by PCS in 
single crystals (from Ref. 8). Lines: Ai {iujn=o) calculated for 
the a and tt bands at T = Tc/4: by solving the imaginary-axis 
Eliashberg equations in the following cases: i) thin solid line: 
pCTir_Q^ and fj, varies with x as shown in Fig.|SJa), open circles; 
ii) dashed line: /x(a;) = /i(0) and F'^'^ varies with x as shown 
in Fig. |SJb) , open squares; iii) thick solid line: both fi and 
F'^'^ vary with x as shown in Fig. |3a) and (b), respectively 
(solid symbols). 



Let us focus for the time being on single crystals. The 
Tc{x) curve (solid circles in Fig. |3J) can be exactly repro- 
duced by adjusting only one of the two free parameters 
of the model, or both of them at the same time (but, in 
this case, the choice of their values is not univocal unless 
one adds another constraint). 

For example, one can keep fi equal to its value in pure 
MgB2 (i.e., ^{x) = /i(0)), and view the substituted com- 
pound as a "disordered" version of MgB2, as if the only 
effect of C substitution was an increase in interband scat- 
tering. This implies neglecting also the phonon harden- 
ing and the electron-doping effects (that actually play 
a leading role in determining the observed Tc{x) curve 
[2^ ) so that the Tc{x) curve is reproduced by only vary- 
ing T""^ . The resulting trend of the interband scattering 
rate is shown in FiglH^a) (open squares). Notice that 
with this approach one cannot obtain critical tempera- 
tures lower than Tc = 25.8 K, that corresponds to the 
isotropic "dirty" limit in which the two gaps merge into 
one of amplitude A = 4.1 meV This is clearly seen in 
the X dependence of the gaps calculated with these val- 
ues of F'^'^, which is reported in Fig. 0] as a dashed line. 
In spite of a rather good agreement between experimen- 
tal and theoretical values of A^, the model predicts an 
increase in A^ which is not observed, and the merging of 
the two gaps at a much lower C content with respect to 
the actual one. 

The opposite case consists in taking into account all 
the effects of substitutions (i.e. phonon hardening and 
electron doping), with no increase of interband scatter- 
ing. In this case, one can keep P'^'^ = 0, and vary /i 
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FIG. 5: (a) Open squares: the x dependence of V'^ necessary 
to reproduce the Tc of single crystals, if iJ.{x) = fi{0). Solid 
squares: the r'^'^{x) curve that allows fitting both Tc and the 
gaps when also fi is varied with x. Lines are only guides to the 
eye. (b) Open circles: the x dependence of the prefactor in 
the Coulomb pseudopotential, /i, that gives the Tc{x) curve 
measured in single crystals (solid symbols in FigEl, when 
po-ir _ g Q^QYid circles: the i-i{x) curve that allows fitting both 
Tc and the gaps (Ao- and A,r) in single crystals, when also 
r"'' is varied. 



with X so as to reproduce the experimental Tc{x) curve. 
With the resulting values of fJ,{x), shown in Figl^^b) as 
open circles, one obtains the x dependence of the gaps 
indicated in Fig. ^ as thin solid lines. It is clear that the 
fj.{x) curve that reproduces the experimental Tc for any C 
content gives values of the large gap that agree rather 
well with the experimental ones (open circles in Fig. ^ 
but gives rise to a decrease in the small gap which is not 
observed experimentally. 

The analysis of the previous two cases suggests that 
the experimental Aa-(x) and AT^ix) curves could be ex- 
plained as due to the interplay between a decrease in /i 
(that makes A^r decrease) and an increase in F"''^ (that 
instead makes Ajr increase). This result has been re- 
cently anticipated by an analysis of the effects of band 
filling and interband scattering j23|. Hence, we will now 
try to fit the experimental x dependence of Tc and of 
the gaps Att and A^r by varying both ^ and T'^'" . The 
best-fitting curves for the gaps (actually, the values of 
A{iuJn=o) at T = Tc/4, for the two bands) are reported 
as thick solid lines in Fig. 0] The choice of the parame- 
ters is univocal, and the resulting x dependencies of T'^'^ 



and fj, are reported as solid symbols in Fig. Ufa) and (b) , 
respectively. 

As shown in Fig.|5ja), the interband scattering remains 
smaller than 2 meV (which is a value reasonable even for 
pure MgB2) up to cc = 0.10 and then increases to make 
the gaps approach each other until they become indistin- 
guishable. The point a.t x = 0.132 in Fig.[3Ja) represents 
the minimum value of T"'^ that gives gap values differ- 
ing less than 0.3 meV (which is approximately the best 
experimental resolution of PCS at 4.2 K). Greater values 
of F'^'^ are allowed as well, since they would give rise to 
gaps even closer to each other. Although the point at 
X = 0.132 might depend on the approximations we are 
using in the present paper, there is no doubt that F'^'^ has 
to increase to reproduce the experimental gap values . 
This increase is thus a general prediction of the two-band 
Ehashberg theory, but its origin in C-substituted MgB2 
is still debated at the moment. According to Ref. |2^ 
carbon substitutions should not change the local lattice 
point symmetry and therefore the interband scattering 
should remain very small as in pure MgB2 pol |. However, 
a (T-vr hybridization might also arise, above x = 0.10, 
from the presence of superstructures or even short-range 
order in the substituted compound 0- It must be said, 
however, that high-resolution TEM has shown no super- 
structures in these single crystals even if the possi- 
bility of short-range order is not ruled out. An alterna- 
tive explanation is based on the observed increase in flux 
pinning and in the normalized resistance on increasing x 
[21j . These effects suggest the existence of microscopic 
defects in the single crystals, acting as scattering cen- 
ters. As indicated by magnetization data, these defects 
might be local inhomogeneities in the C distribution on 
a length scale comparable to ^, that may well give rise 
also to a — TT scattering. 

The values of the Coulomb pseudopotential prefactor, 
/i, that allow reproducing both the Tc and the gap ampli- 
tudes, are reported in Fig. E^b) (solid circles) as a func- 
tion of X. The resulting fi{x) curve is almost linear up to 
X = 0.10, where a change in slope reflects the analogous 
feature of the experimental Tc (see Fig. O. Fig. re- 
ports the values of the components of the Coulomb pseu- 
dopotential matrix, , calculated from eq. [S] by using 
the densities of states {N^{Ef,x) and N^{Ef,x)) from 
density-functional methods, and the values of ^(x) that 
allow best-fitting the experimental gaps (solid symbols in 
Fig. Ob)). It is clear that all the components of the /i* 
matrix monotonically decrease on increasing the C con- 
tent. The large decrease (by a factor of two) of /i or, 
similarly, of fi*^ between x=0 and x=0.1, suggests large 
changes in the electronic screening, that seem to be in- 
compatible with the much smaller changes in the partial 
densities of states (FigQ]). Giving an explanation of this 
puzzle within the two-band model is a very difficult task. 
However, a tentative and qualitative explanation can be 
given in the much simpler single-band case. Let us there- 
fore consider the cr-band quantities alone. Let fi* = fj,*^ 
be the renormalized Coulomb pseudopotential, given by 
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^* = ^[1 + ^ln(£'F/wiog)] ^- Starting from the value 
fi*{x=0) ~ 0.17 (see Fig^ and using Ef=500 meV ^ 
and LU\og—uj-E2g, the value of the bare Coulomb pseudopo- 
tential /i — 0.26 is obtained. From the Morel- Anderson 
definition |2^] of /i, i.e.: 



1 



• In 



f2kF 



(7) 



where ks is the screening wavevector, and using as a 
first approximation the free-electron relationship between 
kp and Ep, one gets fcs(a; = 0)=0.47 A~^. The same 
calculation gives, for x = 0.1, fcs(a; = 0.1)=0.16 A~^, 
so that [{ks{x = 0)/ksix^0.1)f = 8.56. Since in the 
Morel-Anderson model k (where kxp is the 

Thomas-Fermi screening wavevector) and k^p is propor- 
tional to N{Ep), this value has to be compared to the 
ratio N^{Ef,x=0)/N^{Ef,x=0.1)=1.11. The compar- 
ison confirms that the change in the DOS alone cannot 
explain the observed reduction in fi*. However, a large 
increase in the residual resistivity is observed on increas- 
ing the C content [U, so that po{x — 0.l) « 5po{x = 0). 
This suggests that, for some x > 0.1, a metal-to-insulator 
(MIT) transition might be expected. In the hypothe- 
sis that at X = 0.1 the system already lies somewhere 
between the Fermi liquid and the critical regime where 
the MIT occurs, a generalization of the Morel-Anderson 
model has to be used to describe the a; = 0.1 case. 
Within this scenario, fc|(a;=0.1) oc k^p[l + {a/ar)^]~^, 
where r — [pq{x—0.1) — pd/pc, Pc is the critical value of 
the residual resistivity and a, a are constants defined in 
Ref. [2^ Hence one gets 



ksjx = 0) 
ks{x = 0.1) 



= 0) 
N^{x - 0.1) 



l+( — 

ar 



(8) 



from which (a/ar) — 2.78. The parame- 
ter ar expresses the distance from criticality 
(i.e. from the MIT) and can be obtained from 
N^{Ef,x=0.1) = A^^(£:F,a;=0)[l - exp(-ar)], that 
gives ar = 2.3. According to Ref. *2u, this value is 
perfectly compatible with a strongly disordered Fermi 
liquid. Finally, the value of the constant a turns out to 
be a = 6.4 that falls in the range of values allowed in 
Refl2^ and is correctly of the order of the cell parameter. 
In conclusion, the observed drop of /z*^ is due to a 
change in the screening length that, in turns, can be 
justified by the transition to a disordered Fermi liquid 
on increasing the C content. Incidentally, this result 
might further justify the observed increase in interband 
scattering F"''^ at high doping levels. 

At this point, all the parameters entering the two-band 
model in Eliashberg formulation have been determined 
as a function of the C content, so that in principle any 
relevant physical property of the superconducting state 
of Mg(Bi_xCx)2 single crystals can be calculated. For 
the time being, we can calculate the temperature de- 
pendence of the gaps at different C contents, that can 
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FIG. 6: The elements of the Coulomb pseudopotential ma- 
trix, calculated from eq^lby using the densities of states 
from first-principles calculations and the prefactor fi{x) that 
best fits the experimental data (Tc and gaps) in Mg(Bi_xCx)2 
single crystals. 
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FIG. 7: The temperature dependence of ATr{iuj„^o) and 
Acr{iiOri=o) calculated by solving the Eliashberg equations in 
four different cases: x = (solid lines), x = 0.055 (dashed 
lines), 2; = 0.1 (dotted lines) and x — 0.132 (dashed-dotted 
lines) . 



be easily tested by performing PCS measurements as a 
function of temperature. Fig. |7| reports the calculated 
values of A7r(ia;„=o) f^nd Ao-(iw„=o) as a function of T 
given by the solution of the Eliashberg equations in four 
different cases: a; = (solid lines), x = 0.055 (dashed 
lines), X — 0.1 (dotted lines) and x = 0.132 (dash-dotted 
lines). It is worthwhile to notice that, even at high C 
contents, the A^(T) curve shows a negative curvature in 
the whole temperature range as in pure MgB2. This is 
due to the fact that, as shown in Fig. |3 the interband 
coupling does not decrease sensibly on increasing x - oth- 
erwise a positive curvature would be observed in At^{T) 
in the proximity of Tc |27j . 
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FIG. 8: Experimental values of Ao- (open circles) and 
(solid circles) measured by PCS at T = 4.2 K (from Ref. (J 
compared to the values of Ai(iti;„=o) of the a and tt bands 
(lines) calculated by solving the imaginary- axis Eliashberg 
equations at T = Tc/i, when all the physical parameters vary 
with X apart from the interband scattering which is kept equal 
to zero. 

Let us now turn our attention to the experimental re- 
sults obtained in Mg(Bi_xCx)2 polycrystals 7]. As we 
did in the case of single crystals, we start by trying to 
reproduce the experimental Tc{x) curve (open circles in 
Fig. 13) keeping V^"" = and varying the prefactor in the 
Coulomb pseudopotential, /i. Once determined the val- 
ues of fi that give exactly the experimental Tc, we can 
calculate the gaps A(ia;„=o) at T = Tc/4 for the a and 
TT bands. The results are reported as a function of x 
in Fig. IHl (solid lines). Surprisingly, the calculated gaps 
agree very well with those measured by PCS (symbols), 
with no need of interband scattering. This result indi- 
cates that the strong difference between the trend of the 
gaps measured in single crystals ^ and polycrystals W\ 
is very likely to be due to the different nature of the sam- 
ples. Unfortunately, a more detailed discussion would re- 
quire a deeper knowledge of the mechanisms that give rise 
to interband scattering in C-substituted samples, which 
is lacking at the present moment - even though some hy- 
potheses for the increase in T"'^ in single crystals have 



been presented above. 

In conclusion, we have studied the Mg(Bi_xCx)2 sys- 
tem within the effective two-band Eliashberg model, that 
was already shown to be well suited for the description of 
unsubstituted MgB2. In the analysis of the C-substituted 
system, we have used as input parameters the frequency 
of the B-B stretching mode and the partial densities of 
state at the Fermi level, calculated as a function of x 
by first-principles density-functional methods. Adopting 
some reasonable approximations, we have come to a sim- 
plified model with only two adjustable parameters (the 
prefactor in the Coulomb pseudopotential and the inter- 
band scattering rate), whose dependence on x has been 
determined so as to reproduce the experimental values of 
Tc and of the gaps and A^. 

The success of the model in describing the experimen- 
tal findings shows that C-substituted MgB2 is a weak- 
coupling two-band system as the parent compound. In 
details, the results indicate that: i) the experimental be- 
haviour of Tc on increasing x is the results of the de- 
crease in the a — a intraband coupling (due to the filling 
of the a bands) and of a decrease in all the elements of 
the Coulomb pseudopotential matrix, in particular /i*^. 
Without the contribution from fi* , the Tc(x) curve would 
be steeper than experimentally observed |22| : ii) the dif- 
ferent trend of the gaps observed experimentally in sin- 
gle crystals (where the gaps become indistinguishable at 
X — 0.132) and polycrystals (where there is no tendency 
to gap merging) only arises from the different amount of 
interband scattering in the two cases. The increase in 
T°"^ above x ~ 0.10 might arise from short-range order 
in the single crystal structures Q , or from local inhomo- 
geneities in the C content on a microscopic scale |2l| . 

Finally, these results give an indication of what an ideal 
substitution, capable of increasing the Tc of the MgB2 
system, should do, i.e. increase Ao-o-, decrease /i*^., and 
keep the interband scattering as small as in pure MgB2. 
According to eqs. |3| and [S] this is possible if N^{Ep) 
increases and N^{Ep) decreases. 

Many thanks are due to S. Massidda and A. Bian- 
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RBAU01FZ2P and the INTAS Project n.01-0617. 
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